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Preface 

I 

This  study  arose  frcm  the  Air  Force's  need  for  high  efficiency 
solar  cells  for  space  applications.  The  project  was  proposed  by 
Ear  William  P.  Rahilly  of  the  Aero  Propulsion  Laboratory  and  thanks 
are  due  him  for  his  guidance  and  significant  contributions. 
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Abstract 


A  solar  cell  with  a  high  band  gap  window,  a  horo junction,  and  a 
back-barrier  on  a  base  has  been  analyzed  theoretically  under  fairly 
general  conditions.  The  carrier  transport  equations  were  solved  to 
obtain  the  short  circuit  current  density.  The  theoretical  results  are 
applicable  to  similarly  structured  cells  itade  from  other  materials. 
Although  there  is  a  lack  of  data  on  the  properties  of  AlGaAs,  many  of 
the  properties  of  AlGaAs  have  been  modeled  in  this  report  by  drawing 
an  the  work  of  others.  The  AlGaAs  solar  cell  was  modeled  on  a  computer 
and  certain  design  parameters  (such  as  A1  content)  varied  in  order  to 


examine  their  effects 


I  .  Introduction 


Recent  research  on  AlGaAs  (Aluminum  Gallium  Arsenide)  indicates 
that  it  may  be  a  superior  material  for  use  in  solar  cells  (Ref  2; 4) 
AlGaAs  ia^similaf' to  GaAs  in  structure  and  properties.  It  is  antici¬ 
pated  that  fabrication  techniques  already  developed  for  GaAs  cells 
could  be  adapted  to  AlGaAs.  In  this  report  a  particular  multilayer 
AlGaAs  solar  cell  is  analyzed  theoretically,  and  some  of  its  character¬ 
istics  investigated  by  computer  modeling. 

Background 

Single  crystal  GaAs  is  a  III-V  compound  that  forms  in  the  zinc- 
blende  structure.  Ideally  each  Ga  atom  is  surrounded  by  four  As 
atoms,  and  each  As  atom  is  surrounded  by  four  Ga  atoms.  In  AlGaAs 
some  of  the  Ga  atoms  are  replaced  by  A1  atoms.  AlGaAs  is  also  a  III-V 
compound  since  A1  and  Ga  are  both  Column  m  elements  .  The  amount  of 
A1  is  indicated  by  writing  Al-j^Ga^^As,  where  X  is  the  fraction  of  Ga 

sites  of  the  lattice  that  are  occupied  by  A1  atoms. 

\ 

The  primary  effect  of  replacing  some  of  the  Ga  with  A1  is  to 
increase  the  band  gap  energy,  Eg.  As  X  increases  from  0  to  1  Eg 
increases  from  1.424  eV  to  2.168  eV  (Ref  8  ).  As  will  be  shown  later, 
a  larger  Eg  in  an  AIGgAs  solar  cell  leads  to  a  larger  open  circuit 
voltage  but  a  smaller  short  circuit  current.  The  efficiency  of  a 
cell  is  approximately  proportional  to  the  product  of  the  short  circuit 
current  and  open  circuit  voltage.  Calculations  based  on  the  effect 
of  Eg  on  short  circuit  current  and  the  current-voltage  characteristic 
under  AMD  conditions  (earth  orbit  space  sunlight)  at  room  tanperature 
predict  a  maximum  efficiency  for  cells  made  from  materials  with  Eg 


near  1.5  eV(Ref  1:208;3:3  ).  An  optimum  value  of  X  near  .1  was 
expected  for  an  AlxGaj_xAs  solar  cell,  sino3  this  corresponds  to  a 
band  gap,  Eg,  of  about  1.55  eV  (Ref  8).  However,  many  other 
material  parameters  depend  on  X  and  these  also  affect  efficiency.  A 
major  part  of  this  report  is  the  modeling  of  these  parameters. 

Solar  Cell  Structure 

The  structure  of  the  solar  cell  analyzed  and  modeled  in  this  report 
is  shown  in  Figure  la.  The  energy  diagram  is  illustrated  in  Figure  lb. 
There  is  an  anti-reflection  coating  on  the  front.  Regions  1  and  2  are 
p-type  while  regions  3,  4,  and  5  are  n-type,  so  the  cell  is  a  p  on  n 
solar  cell  with  the  p-n  junction  at  Y2-  The  first  layer,  region  1, 
is  a  window.  It  is  thin  and  has  a  high  A1  content  (X  >  .85)  and, 
therefore,  a  high  indirect  band  gap,  such  that  it  will  pass  almost  all 
incident  light.  Its  primary  purpose  is  to  effectively  eliminate  the 
loss  of  photogenerated  minority  carriers  due  to  surface  recombination 
at  the  front  surface.  Even  when  chemically  treated  (application  of  an 
anti-reflective  coating)  there  may  remain  a  large  concentration  of 
surface  states  on  the  front  of  the  cell.  However,  the  high  band  gap 
window  acts  as  a  barrier  to  keep  photogenerated  electrons  in  region  2 
from  reaching  the  front  surface  and  recombining  with  holes.  The  window 
constrains  the  electrons  to  region  2  and  increases  the  probability  of 
their  reaching  the  p-n  junction  at  y'2-  The  width  of  region  1  is 
selected  to  be  . 3u,  large  enough  so  that  tunneling  of  electrons  frcm 
region  2  to  the  front  surface  can  be  neglected,  but  not  too  large 
because  it  is  preferable  to  pass  most  of  the  incident  light  to  regions 
that  are  closer  to  the  junction.  High  band  gap  AlGaAs  windows  are 
already  being  used  on  GaAs  cells  (Ref  15) . 


2 


Regions  2  and  3  are  lower  band  gap  layers  (X  <  .3).  Both  layers 
have  the  same  A1  content  so  that  the  junction  between  them  is  a  homo- 
junction.  Most  of  the  absorption  of  light  and  resultant  photogenera¬ 
tion  of  electron-hole  pairs  takes  place  in  these  two  regions.  The  A1 
content,  X,  is  chosen  to  be  less  than  .3.  This  simplifies  analysis 
because  at  higher  A1  content  (.35  <  X  <  .6)  the  direct  and  indirect 
bands  are  close  together  (Ref  4:CS-5).  The  maximum  efficiency  was 
expected  for  X  near  .1  as  mentioned  previously.  Region  4  has  X  =  .3 
so  that  its  band  gap  is  larger  than  that  of  region  3.  It  will  constrain 
holes  in  region  3  to  move  toward  the  junction.  Region  5  is  the  GaAs 
base  upon  which  the  other  layers  are  grown  and  is  used  because  of 
lattice  matching  and  its  availability  in  single  crystal  form.  Region 
5  has  a  lower  band  gap  than  region  4.  It  will  be  shown  that  region  4 
causes  region  5  to  be  isolated  from  the  other  regions,  so  that  it  has 
no  affect  on  the  characteristics  of  the  solar  cell. 

Organization  of  Report 

The  second  chapter  of  this  report  is  a  fairly  general  derivation 
of  the  short  circuit  current  density  in  terms  of  a  few  parameters  for 
a  solar  cell  which  has  an  energy  diagram  like  that  of  Figure  lb.  Few 
assumptions  are  made  and  the  results  are  applicable  to  similarly 
structured  cells  made  from  other  materials.  In  Chapter  3  the  diode  cur¬ 
rent  is  introduced,  and  the  total  current  as  a  function  of  voltage  and 
the  efficiency  are  examined.  The  treatment  is  still  fairly  general. 

In  Chapter  4  the  parameters  affecting  the  current  and  efficiency 
of  the  solar  cell  are  modeled.  Certain  design  parameters  are  specified 
while  others  are  left  as  variables.  The  model  of  the  solar  cell  is 
then  much  more  specific.  In  Chapter  5  a  brief  description  is  given  of 
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the  program  or  computer  model  of  the  AlGaAs  solar  cell.  The  program  is 
merely  the  compilation  of  the  results  of  Chapters  2,  3,  and  4.  It 
generates  some  of  the  characteristics,  such  as  efficiency,  of  the 
AlGaAs  solar.  Its  input  variables  are  the  design  parameters  which  have 
not  yet  been  specified,  such  as  the  width  of  various  regions  and  the 
A1  content  of  region  2  and  3.  Chapter  5  also  gives  the  results  of 
varying  some  of  these  parameters,  and  these  results  are  examined. 

In  Chapter  6  the  overall  conclusions  and  recorrmendations  are  made. 
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II  Short  Circuit  Current  Density 

T 

The  short  circuit  or  photogenerated  current  density,  J  s>c.,  is 
due  to  the  diffusion  of  charge  carriers,  holes  and  electrons.  The 
carrier  transport  equations,  which  govern  the  net  flow  of  electrons  and 
holes  in  semiconductors,  will  be  derived  in  the  first  section  of  this 
chapter.  In  the  second  section  these  equations  will  be  applied  to  a 
solar  cell  with  the  energy  diagram  given  in  Figure  lb.  The  equations 
will  be  solved  for  J^s.c.,*  In  "the  third  section  the  results  will  be 
analyzed . 

Carrier  Transport  Equation 

In  a  non-magnetic  media  two  of  Maxwell's  equations  may  be  written  as 


V  •  D  =  o  (1) 

V  X  H  =  J  +  |2  (2) 

where 

D  =  electric  displacement  V'3ctor 
H  =  magnetic  intensity  vector 
J  -  current  density  vector 
p  =  free  charge  density 

Taking  the  divergence  of  the  Eqn  (2)  gives 

°  =  v.J  +  -|-7.D  (3) 


since  the  divergence  of  the  curl  of  any  vector  is  zero.  Substituting 
Eqn  (1)  into  Eqn  (3)  yields 

5  •  J  -  -  If  <4> 
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This  is  the  continuity  equation  for  free  charges.  The  volume  integral 
of  v.J  is  the  rate  of  charge  leaving  the  volume,  while  the  volume 

.  ^p 

integral  of  -  —  represents  the  rate  of  decrease  of  charge  in  the 
volume. 

There  are  two  kinds  of  free  carriers  in  the  semiconductor. 
Electrons  move  about  in  the  conduction  band  and  holes  are  transported 
in  the  valence  band.  Thus 


P  =  q  (P-N)  (5) 

where 

q  =  the  elementary  charge  (1.6  x  10“^  coulombs) 
p  =  density  of  holes  (g/cmP) 

N  =  density  of  electrons  (#/cm3) 


The  current  density,  J,  can  be  broken  up  into  terms  for  the  drift  and 
diffusion  of  holes  and  elections 


(6) 


The  drift  current  densities  come  from  Ohm's  Law  and  are  due  to  the 
net  drifting  of  charge  carriers  in  an  electric  field  E. 


J§rlft  "  <S‘'pEP  <7> 

Jdrift  =  cll,nEN 

where 

Pp  =  mobility  of  holes  (cm^/volt-sec) 
pn  =  mobility  of  electrons  (cm^/volt-sec) 

The  diffusion  current  densities  are  based  on  Fick’s  Law  and  are  due  to 
the  average  flow  from  a  region  of  high  concentration  to  region  of 
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low  concentration.  In  symbolic  form 


JP  = 
diff 

-  qDp  VP 

(9) 

tN 

Jdiff 

qDn  VN 

(10) 

where 

Dp  =  diffusivity  of  holes  (cm^/sec) 

Dn  =  diffusivity  of  electrons  (cm^/sec) 

Eqn  (4)  may  now  be  rewritten  using  Egns  (5)  through  (10)  as 

V  •  (ppEP  -  Dp  VP  +  UjjEN  +  Dn  VN)  =  ||  -  J|  (11) 

where  q  has  been  divided  out.  In  order  to  obtain  separate  equations 
governing  the  transport  of  electrons  and  holes  from  Eqn  (11)  carrier 
generation  and  recombination  terms  must  be  included.  This  is  because 
in  general  the  generation  and  recombination  rate  of  holes  and  electrons 
are  a  function  of  both  P  and  N.  The  Shockley-Hall  carrier  generation- 
recombination  rate  is  given  by 

U  =  Z  ( i  o\ 

(^Pt)Tno  +  (WO'po 

where 

n^  =  intrinsic  carrier  density  (cnf^) 

Pt,Nt  =  trap  densities  (cnT^) 

Tno.Tpo  =  carrier  lifetimes  (sec) 

as  defined  by  Sze  (Ref  12:48).  Note  that  since  holes  and  electrons 
are  oppositely  charged  and  since  both  generation  and  recombination 
involve  electron-hole  pairs,  neither  generation  nor  recombination 
affect  the  total  free  charge  density,  p.  Since  the  carriers  move 
independently  in  their  respective  gands,  Eqn  (11)  may  be  rewritten  as 
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se  i  m 


t 


» 


|f£  +  V- (ijpEP  -  DpVP)-G  +  U^j 


||  +  V.(pn  EN  +  Dn  VN)  +  G 


-  U  J  =  0 


(13) 


and  each  straight  bracketed  term  can  be  set  equal  to  zero.  In  the 
present  application  the  generation  rate,  G,  is  taken  to  be  a  function 
of  position  only.  Thus  Eqn  (13)  fragments  into 


+  V*(yp  EP  -  Dp  VP)  -  G  +  U  =  0  (14) 

3N  +  V-(-pn  EN  -  Dn  VN)  -  G  +  U  =  0  (15) 

at 

E ins  (14)  and  (15)  are  the  carrier  transport  equations  for  holes  and 
electrons  respectively.  They  are  coupled  by  the  recombination  rate, 

U,  given  in  Eqn  (12)  and  the  electric  field  vector  E.  Using  the 
relation  D  =  eE,  where  e  =  permittivity,  in  Eqn  (1)  gives  the  third 
equation  needed  to  determine  E,  P,  and  N.  However,  the  electric  field 
will  be  assumed  to  be  zero  so  that  only  Eqns  (14)  and  (15)  are  needed. 


Solution  of  the  Transport  Equations 

If  the  solar  cell  is  unilluminated  it  is  merely  a  diode.  In  a 
given  region  and  with  no  external  voltage  bias,  the  electron  and  hole 
concentrations  will  be  at  their  thermal  equilibrium  values  No  and  PD 
respectively.  The  short  circuit  current  density  will  be  zero.  If  the 
cell  is  illuminated  the  concentrations  are  increased  to  N  =  No  +  n 
and  P  =  Po  +  p,  where  n  and  p  are  the  excess  electron  and  hole  con¬ 
centrations  due  to  photogeneration.  The  short  circuit  current  or 
light  generated  current  is  due  to  the  excess  carriers.  In  order  to 
obtain  Js.c.  the  excess  minority  carrier  concentrations  in  each  region 
must  be  found. 
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In  Figure  lb  the  energy  band  edges  for  the  electrons  and  holes 


change  near  the  edges  of  each  region  but  are  constant  in  between. 

The  energy  bands  will  be  assumed  to  change  abruptly  at  each  interface, 
except  in  Chapter  3  where  the  carrier  depletion  width  at  the  p-n 
junction  will  be  used  to  obtain  the  diode  current.  Actually  the 
energy  bands  change  over  a  snail  width  and  an  electric  field  (pro¬ 
portional  to  the  gradient  of  the  energy  band  edge)  exists  there. 

Since  the  energy  bands  are  constant  with  position  in  each  region  no 
electric  field  exists.  Therefore,  in  a  given  region  assuming  time 
independence  Eqn  (15)  can  be  written  as 

DnV2N  +  G-  U  =  0  (16) 


Photogeneration  is  the  interaction  of  a  photon  and  an  electron 
in  the  valence  band  which  raises  the  electron  to  the  conduction  band 
and  leaves  a  hole  in  the  valence  band.  It  is  assumed  that  no  other 


absorption  processes  (such  as  free  carrier  absorption)  are  insignifi¬ 
cant  for  the  wavelengths  of  interest  (.35u  to  .90p).  At  a  particular 
wavelength  of  light  the  generation  rate  G  is  proportional  to  the  flux 
of  photons  at  that  wavelength.  However,  the  photon  flux  decreases 
with  increasing  y  due  to  absorption. 


Let  Fj_  be  the  flux  of  photons  at  X  incident  on  the  very  front  of 
region  1.  The  photon  flux  in  region  1  can  be  written  as  F^e-01!^ , 


where  is  the  absorption  coefficient  at  X.  Since  the  generation  rate 


of  carriers  is  equal  to  the  absorption  rate  of  photons,  from  Beer's  Law 


Gi  =  «x  Fx  e-^ 


(17) 
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The  y  dimension  of  each  layer  of  the  solar  cell  is  much  smaller  than 
its  other  two  dimensions,  and  since  the  generation  rate  varies  only  with 
y,  the  carrier  densities  wall  be  assumed  to  vary  only  with  y.  In 
addition,  since  NQ  is  constant  in  a  given  region  the  variation  of  N 
with  y  is  due  solely  to  n  the  excess  electron  concentration.  In 
region  1  Eqn  (16)  can  now  be  written  as 

°ni  f~J  -  UX  =  -  aiFie-aiy  (17) 

In  a  p-type  region  P0»No.  The  excess  carriers  due  to  photo¬ 
generation  have  little  effect  on  P  compared  to  the  effect  on  N. 

P  =  P0+p«P0  (18) 

N  =  Nq  +  n  (19) 

Since  the  term  Prno  dominates  in  the  denominator  of  Eqn  (12)  and  since 
nf  =  PqN0,  the  recombination  rate  in  region  1  reduces  to 

Ui^2I  (20) 

tnoi 

Substituting  into  Eqn  (17)  yields 

Dbopi  -  a  ,  -  aiFie-0iy  (2i> 

The  diffusion  length  of  a  carrier  is  related  to  the  diffusivity  and 
lifetime  by 

L  =  Dt  (22) 

Therefore  Eqn  (21)  can  be  written  as 

“"1^  -  =  -ajFje-n*  (23) 
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Equations  similar  to  this  may  be  written  for  regions  2,  3,  and  4. 


°n2  =  -  ct2F2e-a2(y-yi)  (24) 

n2 

Dps  IP  -  Sgp  -  -  «3F3e-“3(y-y2)  (25) 

°P4  IP  -  ^  *  -  “A8'*41™*  (26) 

where  =  Fj^  e-ak-l(yk-l”yk-2)  is  the  photon  flux  incident  on  the  front 
of  region  k.  Eqns  (23)  through  (24)  along  with  boundary  conditions  not 
yet  introduced  determine  the  excess  minority  carrier  concentrations 
in  regions  1  through  4.  The  term  on  the  right  hand  side  of  each 
equation  is  the  forcing  term.  It  determines  the  particular  solution 
to  the  differential  equation. 

Since  it  is  the  minority  carriers  in  each  region  that  are  of 
primary  interest,  the  subscripts  p  and  n  will  be  dropped  when  only 
minority  carriers  are  being  discussed. 

The  solutions  to  Eqns  (23)  through  (26)  are 

ni(y)  =  Ai  cosh  (^£^0  +  B}  sinh  +  C2e-Cliy  (27) 

n2(y)  =  A2  cosh  (^p)  +  B2  sinh  C^2)  + 

P3(y)  =  A3  cosh  (^p)  +  B3  sinh  (y~|2)  +  C3e_°l2(y-^) 

P4(y)  =  A4  cosh  (^qp2)  +  B4  sinh  (^qp)  +  c4e_a^^y_f3o) 


where 
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The  rate  of  excess  electrons  coming  out  of  the  surface  on  the  right  side 
must  equal  the  rate  going  in  on  the  left  side  minus  the  rate  at  which 
they  are  lost  due  to  interface  recorrbination  states.  Therefore 

Da  3s|ai  .  ^  dalai  +  Sl  „2(n)  (34) 

At  Y2,  the  junction,  another  electric  field  exists  which  ties  the 
excess  electron  concentration  to  zero  on  the  left  side. 

n2(y2>  =  0  (35) 

Using  similar  arguments  similar  B.C.s  result  for  the  excess  hole  density 
on  the  p  side  of  the  solar  cell; 


P3(y2)  =  0 

(36) 

n  =  03  itep)  +  S3  P3(y3) 

(37) 

P4(y3>  =  0 

(38) 

P4(y4)  =  0 

(39) 

Hie  integration  constants  for  ?4(y)  are  determined  by  Eqns  (38)  and  (39). 
Hie  hole  density  in  region  4  is  affected  by  the  other  regions  only  through 
which  is  affected  by  the  absorption  of  light  in  the  preceding  regions. 
Pq(y)  affects  P3(y)  through  Eqn  (37).  Similarly  nj(y)  is  independent 
of  the  other  regions  but  it  affects  n2(y)  through  Eqn  (34). 

The  following  definitions  are  made  for  convenience: 

zk  =  yk  -  yk-1  =  width  of  region  k  (40) 

Ux  =  cosh  (-  gjO  ,  V2=  sinh  (-  g)  (41-42) 

U2  =  cosh  (-  ,  V2=  sinh  (-  (43-44) 
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U3  =  cosh  (|&)  ,  V3  =  sinh  (^)  (45^46) 

U4  =  cosh  (jJ)  ,  V4  =  sinh  (|j)  (47^8) 

U  here  should  not  be  confused  with  the  recombination  rate  given  in 
Eqn  (12). 

The  direction  of  positive  current  will  be  taken  as  the  negative  y 
direction.  The  following  current  densities  are  pertinent: 


Ji  =  -q  Dj  =  -Q  Dl  (§-  -  a-^e-'*!2!)  (49) 

J2  =  -<1  d2  ^dy^  =  °2  (^  "  a2C2e-Cl2z2)  (50) 

J3  =  Q  D3  dP3|y-2-)  =  q  D3  (||  -  a3c3)  (51) 

J4  =  q  °4  «  q  Z>4  (§  -  «4 C4)  (52) 


J 1  is  the  electron  current  density  at  yj.  It  represents  the  injection 
of  electrons  from  region  1  into  region  2.  J3  is  the  electron  current 
density  at  y3-  It  represents  the  electrons  flowing  into  the  junction 
from  region  2.  J3  is  the  hole  current  density  at  y3-  It  represents 
the  holes  flowing  into  the  junction.  Finally,  J4  is  hole  current 
density  at  y3.  It  represents  the  injection  of  holes  from  region  4  into 
region  3. 

The  and  coefficients  can  now  be  presented  in  a  way  that  is 
actually  much  less  painful  than  might  otherwise  be  possible.  The  result 
of  applying  the  B.C.s,  Eqns  (32)  through  (34)  and  the  above  definitions 
is: 


A]_  =  -  Ci  e_alzI 


(53) 
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Bl 


=  CKSo+Diai  +  (5gl  _  s0U1)e-°tlzl) 
/DlUl  c  \ 

(TT  "  s°Vl) 


(54) 


A2  =  -  C2  e~a2z2 


B2  =  ~  ~  +  C2  (Si  +  D2a2  +  (^  -  SiU2)e-a2z2) 

P2U2  -  SlV2 

L2 


(55) 

(56) 


A3  =  -C3 

B3  =  ^  +  C3  ((-S3  +  D3  a3)e-a3r3  +  +  S3U3) 


D3U3 

“IT  +  SsV3 


A4  =  -  c4 

B4  =  C4  (U4  -  ea4^) 
V4 


(57) 

(58) 


(59) 


(60) 


These  are  the  coefficients  that  were  needed  in  Eqns  (27)  through  (30). 
The  excess  minority  carrier  concentrations  in  regions  1  through  4  have 
now  been  determined.  However,  it  is  the  current  densities  that  are  of 
primary  interest.  Using  Eqns  (49)  through  (60)  the  current  densities 
may  be  written  more  explicitly  as 


Jl  =  -  q  D1C1 


So*i“1  +  tog1  -  SqUi^i  __  Zl 


Dlul  ~  S0V1L1 


-aje 


(61) 


J2  =  -  qD2 


+  c2  ^1  +  D2^2  +'<¥0^  ~  S1U2)e'a2Z2) 


-a2C2e'^t2z2 


D2U2  -  SiV2L2 


(61) 
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(63) 


J3  =  qD3  ^  +  C3  ((-S3  +  D3«3)e_ct6z3  +  +  s3U3) 

D3U3  +  S3VsL3 


TT/1  —  o —  a42vl 

J4  =  q  D4C4  (■  4  ~  - a4)  (64) 

The  preceding  derivation  was  done  for  only  one  wavelength  of  light, 

A.  Since  the  absorption  coefficients  are  a  function  of  A,  the  spectral 
photon  fluxes  are  needed.  If  F^  =  F]<( A)  is  the  spectral  photon  flux  incident 
on  the  front  of  region  k  then  the  preceding  derivation  gives  the 
spectral  excess  minority  carrier  densities,  nj^( A )  and  pk( A ) ,  and  the 
spectra]  current  densities  J^( A) .  The  total  excess  minority  carrier 
concentrations  and  current  densities  can  be  obtained  by  integra-ing 
over  wavelength. 

T 

nk  =  -(nk(A)  dA  (65) 

Jk  =  /Jk(A)  dA  (66) 

T  T 

J3  is  the  electron  current  density  at  the  p-n  junction.  J3  is  the 

hole  current  density  at  the  p-n  junction.  The  total  current  density 
in  the  junction  is  their  sum.  Since  time  independence  was  assumed  the 
carrier  concentrations  must  remain  constant  with  time.  Therefore  the 
total  current  density  is  the  same  everywhere  in  the  cell;  it  does  not  vary 
with  y.  The  short  circuit  or  light  generated  current  is 

JJ.O.  -  +  J3  <67> 
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Analysis  of  Theoretical  Results 

The  spectral  current  densities  given  in  Eqns  (61)  through  (64) 
are  of  similar  form.  The  term  ‘*'S  conir^u~ 

tion  to  J2  due  to  the  current  injection  from  region  1  into  region  2. 

This  contribution  is  less  than  the  injection  current  density,  Ji ,  since 
U2  >  1  and  V2  <  0  while  D2.  Sy,  and  L2  are  positive.  The  same  is  true 
for  the  contribution  of  J4  to  J3.  This  is  because  some  of  the  injected 
carriers  are  lost  through  recombination  either  at  an  interface  or  in  the 
bulk  of  the  region.  The  bulk  recombination  is  represented  by  the 
diffusion  lengths  L^. 

Except  for  the  contribution  due  to  current  injection,  J2  has 
exactly  the  same  form  as  Ji.  J2  and  J3  are  very  similar.  The  difference 
in  the  form  of  Eqns  (62)  and  (63)  are  due  to  the  direction  of  light, 
toward  the  junction  in  region  2  and  away  from  the  junction  in  region  3. 

As  previously  mentioned,  the  excess  hole  density  in  region  4  is 
independent  of  region  5.  Since  it  is  the  slope  of  the  excess  hole 
density  at  yo  that  determines  J4,  jjj  is  independent  of  region  5.  There- 
fore  J3  and  Jg.c<  are  also  independent  of  region  5.  This  is  due  to  the 
B.C.  at  y4,  Eqn  (39),  which  is  a  result  of  holes  in  region  4  having  a 
higher  potential  energy  than  holes  in  region  5.  The  primary  purpose  of 
region  4  is  to  act  as  a  barrier,  preventing  holes  in  region  3  from  moving 
beyond  V3  and  therefore  constraining  them  to  move  toward  the  p-n  junction. 
However,  region  4  has  other  effects.  Another  beneficial  effect  is  the 
injection  of  holes  from  region  4  into  region  3.  The  injection  current 
density,  J^,  will  probably  be  snail  for  a  practical  solar  cell  because 
much  of  the  light  will  be  absorbed  before  it  reaches  region  4.  Still, 
several  steps  can  be  taken  to  optimize  the  current  injection.  J4  can 
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be  increased  by  decreasing  the  AI  content  in  region  4,  which  increases 
the  absorption.  However,  the  band  gap  in  region  4  must  be  greater  than 
the  band  gap  in  region  3  in  order  for  region  4  to  act  as  a  barrier  to 
the  holes  in  region  3.  The  band  gap  increases  with  Al  content.  The 
Al  content  in  region  4  must  be  somewhat  higher  than  the  Al  content  in 
region  3  in  order  for  the  B.C.s  at  y3,  Eqns  (37)  and  (38),  to  hold,  and 
in  order  for  region  4  to  act  as  a  barrier. 

As  mentioned  in  Chapter  1,  it  was  anticipated  that  the  greatest 
efficiencies  would  be  achieved  for  an  Al  content  in  regions  2  and  3  of 
^2,3  near  .1.  The  Al  content  in  region  4  was  chosen  to  be  X4  =  .3. 
Therefore,  the  B.C.s  at  y3,  Eqns  (37)  and  (38),  will  be  good  approxima¬ 
tions  for  X23  up  to  about  .25.  A  higher  X4  would  produce  a  lower  J4. 
Also,  as  indicated  in  Chapter  1,  the  direct  and  indirect  bands  are 
close  together  for  .35  <  x  <  .6  so  that  modeling  of  seme  material 
parameters  becomes  very  difficult  at  those  Al  contents  since  significant 
electron  transport  occurs  in  three  conduction  bands. 

J4  can  also  be  optimized  by  adjusting  the  width  of  region  4,  z^. 
This  could  be  done  by  varying  Z4  in  the  computer  model.  A  large  Z4  is 
probably  desirable,  since  the  electric  field  at  y4  acts  as  a  drain  on 
excess  holes  in  region  4.  This  drain  would  then  be  far  removed  from 
the  front  of  the  region  where  the  photon  flux  is  higher  and  where  most 
of  the  photogeneration  in  the  region  takes  place. 

Region  4  also  has  an  adverse  affect.  This  is  due  to  the  recombina¬ 
tion  at  y3,  which  acts  as  a  drain  on  the  excess  holes  in  region  3.  In 
the  worst  case,  when  S3  -*  ®,  Eqn  (63)  reduces  to 

U3-e-a3z3 

J3  =  q  D3C3  (  -3-g-T-- - ct3)  (68) 
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If  region  4  were  not  included  in  the  solar  cell  then  regions  3  and 


5  would  be  adjacent  to  each  other.  Region  3  has  a  higher  A1  content, 
larger  band  gap,  and  higher  hole  potential  energy  than  region  5.  There¬ 
fore  the  electric  field  at  their  interface  would  be  in  the  positive  y 
direction.  One  of  the  boundary  conditions  tbre  would  be 

P3(Y3)  =  0  (69) 

Thus  the  excess  liole  density  in  region  3  would  be  tied  to  zero  at  both 
edges.  This  was  the  case  for  P4(y)  before.  J3  would  then  have  the  same 
fonn  that  J4  had  in  Eqn  (64).  But  J3  in  Eqn  (68)  is  also  of  that  same 
form.  J3  is  the  same  for  both  cases,  changing  the  B.C.  at  y3  to 
Eqn  (69)  or  letting  S3  become  infinite.  Therefore,  J3  and  Jw  r  will  be 
at  least  as  large  with  region  4  included  in  the  solar  cell  as  without 
it. 


* 
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Ill  Diode  Current  Density  and  Efficiency 


In  this  chapter  the  diode  current  density,  Jg,  is  introduced.  The 
total  current  density,  jJ  and  the  efficiency  of  the  solar  cell  are  also 
discussed . 

The  depletion  width,  W,  of  the  junction  was  neglected  in  determining 
Js.c.'  W  is  very  small  and  little  photogeneration  takes  place  in  the 
junction.  However,  W  is  important  in  determining  Jg.  If  dopants  with 
low  diffusion  coefficients  are  used,  and  the  doping  levels  on  both  sides 
of  the  junction  are  high,  the  junction  can  be  assumed  an  abrupt  junction. 
The  depletion  width  is  then  given  by  (Ref  12:89). 

w  -  ¥  >  (VB-V>  <7t» 

where 

eg  =  static  permittivity  (farad/m) 

P20  =  thermal  equilibrium  hole  concentration  in  region  2 
N30  =  thermal  equilibrium  electron  concentration  in  region  2 
V  =  voltage  across  the  solar  cell,  front  to  back 
Vg  =  built-in  voltage  of  the  junction 

P20  is  approximately  equal  to  the  concentration  of  ionized  acceptors 
in  region  2.  N30  is  approximately  equal  to  the  concentration  of 

ionized  donors  in  region  3. 

The  built-in  voltage  is  given  by  (Ref  12:87) 

vB  =  m  (P20N30)  (71) 

i2 


1 


where 
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n 


kg  =  Boltanan's  constant  (1.38  x  10”23  j/°k) 

T  =  absolute  temperature  (°K) 

=  ni3  =  the  intrinsic  carrier  density  of  regions  2  and  3 

The  current  due  to  electron-hole  recombination  in  the  depletion 
region  is  assumed  to  be  the  dominant  dark  current.  According  to  the 
Sah-Noyce-Shockley  theory  (Ref  10),  the  diode  durrent  density  has  a 
maximum  value  of 


JD 


q  ni2  W 
L2L3 


D2D3  2  sinh  (qV/2kBT) 
-  •  q  TVq-V )7kB* 


2 


This  can  be  approximated  by 


(72) 


JD  =  q  ni2  W  D2°3  (e^/^BT.i)  (73) 

2  L2L3 

where  1<  a  <2. 

The  worst  case,  a  =  2,  will  be  assumed. 

The  diode  current  is  defined  as  going  from  the  p  side  to  the 
n  side  in  the  diode  or  cell,  the  positive  y  direction  in  the  cell  of 
Figure  2.  However,  the  photogenerated  or  short  circuit  current  and  the 
total  current  are  defined  as  going  from  the  n  side  to  p  side  within 
the  cell.  Therefore,  the  total  current  density  is  given  by 

jT  =  JL.  "  JD  <74) 

If  the  front  and  back  of  the  solar  cell  are  connected  by  a 

resistanceless  wire  (short-circuited)  then  V  must  be  zero.  Therefore, 

JD  =  0  and  the  total  current  density,  J  ,  is  equal  to  the  photo- 

T 

generated  or  short  circuit  current  density  Jg.c#- 
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If  the  front  and  back  of  the  solar  cell  are  unconnected  (open- 
circuit)  then  jT  =  0,  and 

jJ.C.  -  JD  (75) 

Jg  c.  is  determined  by  the  illumination  and  is  independent  of  voltage. 
Eqn  (75)  determines  the  open -circuit  voltage  Voc>,  since  Jp  is  a 
function  of  voltage. 

Under  normal  operation  a  load  of  resistance  R  is  placed  across  the 
cell;  the  cell  acts  as  the  source.  The  condition 

V  =  1TR  =  JtAR  (76) 

where  A  =  area  of  the  solar  cell 

must  be  satisfied.  Eqn  (76)  determines  the  voltage,  V,  and  the  current 
density,  JT,  of  the  solar  cell  since  ^  is  a  function  of  voltage.  As 
R  goes  from  0  to  «  ^  goes  frcm  Js.c.  to  0  and  V  goes  from  0  to  VQ  c> . 

A  typical  current  density  versus  voltage  curve  for  a  solar  cell  is 
shown  in  Figure  2. 

The  power  output  Pout  of  the  cell  is  given  by 

pout  =  V  (77) 

The  efficiency  n  is  given  by 

n  =  Soul  =  JTV,  .  (78) 

n  Pin  Pin/A 

For  a  given  illumination  Pfjn/A  is  a  constant.  Therefore  the  maximum  n 
and  maximum  Pout  occur  when  the  product  of  and  V  is  a  maximum.  This 
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CURRENT  DENSITY  (mA/cm 


VOLTAGE(volts) 


Figure  2.  Typical  J-V  Characterl stl c  of  a  Solar  Cell 
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optimum  operating  point  is  indicated  by  the  dotted  lines  in  Figure  2. 
n  at  this  point  is  usually  referred  to  as  the  efficiency  of  the  solar 
cell. 

T 

Note  that  decreasing  the  photogenerated  current  density  Js.c.  shifts 
the  curve  in  Figure  2  down. 
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IV  Operating  Conditions,  Design  Parameters, 
and  Material  Character ist ics 

The  preceding  chapters  give  the  gross  characteristics  of  the  solar 
T 

cell,  such  as  Jsc_  and  n ,  in  terms  of  many  parameters  whose  values  have 
not  yet  been  specified.  Some  of  these  parameters  are  operating  conditions 
such  as  temperature.  The  conditions  under  which  the  model  solar  cell 
is  assumed  to  be  operating  will  be  specified  in  this  chapter.  Some  of 
the  parameters  are  design  parameters.  It  is  important  to  note  that 
some  design  decisions  were  already  made.  They  are  implied  by  structure 
and  energy  diagram  of  Figure  1,  by  the  boundary  conditions  used  in 
Chapter  2,  and  by  other  assumptions  already  made.  These  require  certain 
limits  on  the  design  parameters.  These  limits  are  clearly  stated  in 
this  chapter.  Some  design  parameters  are  left  as  variables  that  can  be 
altered  on  different  computer  runs.  Other  design  parameters  will  be 
specified  in  order  to  reduce  the  complexity  of  the  computer  model. 

Many  of  the  material  characteristics  depend  on  design  parameters. 

The  doping  levels  and  A1  contents  of  the  various  regions  are  design 
parameters  because  they  affect  the  material's  electrical  and  light  ab¬ 
sorption  properties.  The  material  properties  of  AlGaAs  used  in  the 
computer  model  are  specified  or  given  in  terms  of  design  parameters. 

This  is  the  most  difficult  and  least  accurate  task  attempted  in  this 
report.  The  great  advantage  of  ternary  compounds  like  AlGaAs  is  that 
they  give  the  designer  the  choice  of  a  spectrum  of  materials  and 
material  characteristics.  However,  this  poses  the  problem  of  determining 
the  characteristics  of  the  material  as  a  function  of  the  added  dimension. 
Modeling  of  9ome  material  properties  is  also  difficult  because  they 
are  fabrication  dependent. 
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The  solar  cell  is  assumed  to  have  a  temperature,  T,  of  300°K. 
Most  information  on  AlGaAs  is  specified  at  this  temperature.  This  is 
also  the  temperature  at  which  results  for  solar  cells  are  ccrrmonly 


reported. 

Incident  Light 

The  light  incident  on  the  cell  is  specified  as  sunlight  in  space 
near  the  earth  at  normal  incidence,  since  the  cell  is  intended  for  space 
applications.  These  conditions  are  referred  to  as  AMO  (Air  Mass  Zero) 
conditions.  The  solar  spectral  ir radiance,  P( X) ,  for  AMO  conditions 
proposed  by  Thekaekara,  et  al  (Ref  14:74)  is  used  in  this  study. 

P(X)  is  shown  in  Figure  3. 

The  total  irradiance  or  power  per  unit  area  was  determined  by 
Thekaekara,  et  al  (Ref  14:75)  to  be 

oo 

Pin  =  /  P(X)  dx  =  .13530  watts/cm2  (79) 

A  o 

This  is  used  in  Eqn  (78)  to  determine  the  efficiency  of  the  solar  cell. 
The  energy  of  a  photon  of  wavelength  x  is 

E  =  ^  (80) 

where  h  =  Planck's  constant  (6.626xl0-2^  joule  sec) 
c  =  speed  of  light  (2.998x10®  m/sec) 

The  spectral  photon  flux  incident  on  the  solar  cell  is  then 

F0(x)  =  (81) 
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Figure  3.  Solar  Spectral  Irradiance  for  Air  Mass  Zero 
Conditions  (Ref  14:74) 


* 
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Reflectivity 


An  uncoated  GaAs  solar  cell  with  an  AlGaAs  window  has  a  reflectivity, 
R(A),  greater  than  .30  at  all  wavelengths  of  interest  (Ref  9  ).  A  two- 
layer  anti-reflection  coating  has  reduced  R(A)  to  between  .02  and  .05 
for  A  between  .4a  and  .88a  (Ref  9).  In  the  present  study,  R(x)  is 
assumed  to  be  .03  for  .35a  <  X  >  .9a  and  1  at  all  other  wavelengths. 

This  light  passed  corresponds  to  photon  energies  of  1.378eV  <  E  < 

3.542  eV.  If  longer  wavelength  photons  were  passed  into  the  cell  they 
would  produce  almost  no  photogenerated  carriers,  since  their  energies 
are  less  than  the  band  gap  energies  of  all  regions.  If  higher 
frequency  photons  were  passed  into  the  cell  the  effect  would  also  be 
anall,  since  photon  flux  for  X  <  .35a  is  small. 

The  spectral  photon  flux  passed  into  region  1  is 

F]_( A)  *  (1  -  R(X»  F0  (X)  (82) 

The  spectral  photon  flux  entering  each  of  the  other  regions  is 
determined  by  Eqn  (31b). 

Widths  of  Regions 

The  width  of  each  region  is  a  design  parameter  and  is  left  as  a 
variable  in  the  computer  model.  However,  there  are  constraints  on  these 
variables.  Since  quantum  mechanical  tunneling  of  minority  carriers 
through  the  barriers  (regions  1  and  4)  was  neglected  in  the  derivation 
of  Jg  c  in  Chapter  2,  the  width  of  these  barriers  must  be  large  enough 
to  make  tunneling  negligible.  In  addition,  practical  limits  exist  on 
how  thin  a  layer  can  be  made.  The  constraints  used  for  the  region 
widths  are  indicated  in  Table  I. 
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Table  I 


Some  of  the  Parameters  used  in  the  Computer  Model 


Parameters  (Units) 

Region 

1 

2 

3 

4 

Z  (p) 

Zi  >  0.3 

Z2  >  0.3 

Z3  >  0.3 

Z4  >  1.0 

X 

.9 

X2  <  .25 

X3  =  X2* 

.3 

Er  (eV) 

2.779 

Bqn  (83) 

* 

1.798 

Ex  (eV) 

2.128 

Eqn  (84) 

* 

1.950 

Er,  (eV) 

3.66 

Eqn  (85) 

* 

3.06 

c<(\)  (cm--*-) 

Eqn  (87) 

Eqn  (87) 

* 

Eqn  (87) 

% 

.79  + 

Eqn  (90)+ 

* 

.092 

Mp 

.65 

.65 

*  + 

.65  + 

n-j_  (cm-3) 

- 

Eqn  (96) 

* 

- 

KS 

11.1 

Eqn  (98) 

* 

12.4 

Kd 

8.7 

Eqn  (99) 

★ 

10.1 

Dopant 

Be 

Be 

Te 

Te 

Ej_  (eV) 

.0424 

Eqn  (102) 

Eqn  (101) 

.00467 

Po  (1017cm-3) 

7.095 

Eqn  (117) 

- 

- 

Nq  (1017cm-3) 

- 

Eqn  (115) 

4.13° 

L  (u) 

2 

2 

2 

2 

p  (cm2 /V- sec) 

99.4 

Eqn  (119) 

Eqn  (121) 

164 

D  (cm^/sec) 

2.57 

Eqn  (122) 

Eqn  (122) 

4.23 

+  =  relative  mass  of  minority  carrier 
*  =  same  value  as  in  region  2 
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If  Z]_  is  increased  above  0.3p  there  is  slightly  more  photogeneration 
in  region  1  and  less  in  the  more  favorable  regions  next  to  the  p-n 
junction,  so  Js.c.  and  n  are  decreased.  Therefore,  Zi  is  held  at  0.3u 
in  all  computer  runs,  although  it  is  left  as  a  variable  in  the  program 
since  it  is  easy  to  do  so.  For  practical  solar  cell  dimensions  little 
photogeneration  will  take  place  in  region  4.  Z4  has  little  affect  on 
cell  efficiency  as  long  as  it  is  large  enough  to  prevent  significant 
tunneling  of  holes  from  region  3  to  region  5.  Z2  and  Z3  are  the 
dimensions  of  primary  interest.  Their  affect  is  examined  in  the  next 
chapter . 

A1  Content  of  Regions 

As  discussed  in  Chapter  2,  the  relative  A1  content  of  the  various 
regions  must  be  chosen  such  that  the  B.C.s  used  in  deriving  the 
photogenerated  current  are  valid.  The  A1  content  of  regions  1  through 
4  are  indicated  in  Table  I.  The  A1  content  of  regions  2  and  3  are  the 
same  (X2  =  X3)  and  are  left  as  a  single  program  variable. 

Energy  Band  Edges 

In  AlGaAs  the  maximum  in  the  two  highest  valence  bands  occurs 
where  the  momentum  vector,  K,  is  zero.  Two  of  the  three  valence  bands 
are  degenerate  at  K  =  0  and  the  split -off  band  is  ignored.  Three 
minima  occur  in  the  lowest  conduction  band.  The  electron  potential 
energy  level  at  these  minima  relative  to  that  of  the  maxima  in  the 
valence  band  are  Ef,  at  K  =  0,  E^,  at  the  minima  along  the  ^100^  lattice 
direction,  and  E[ ,  at  the  minima  along  the  (ill)  lattice  direction. 

These  energies  in  eV  vary  with  A1  content  (Ref  8  ). 
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Ep  =  1.424  +  1.247  X  eV  for  X  <  .45 

1.424  +  1.247  X  +  1.147  (X-0.45)2eV  for  X  >  .45 

(83) 

Ex  =  1.900  +  0.125  X  +  0.143  X2  eV  (84) 

El  =  1.708  +  0.642  X  eV  (85) 

Figure  4  shows  the  dependence  of  the  band  edges  on  A1  content. 

AlGaAs  has  a  direct-indirect  crossover  at  X  =  0.45.  For  X  <  .45 
the  material  is  a  direct  semiconductor,  since  Ep  is  the  lowest  energy. 
For  X  >  .45  the  material  is  an  indirect  semiconductor,  since  Ex  is  the 
lowest  energy.  For  the  A1  contents  of  interest  the  L  minimum  has  little 
affect  and  will  be  ignored  because  it  is  an  indirect  minima  and  because 
El  is  always  sufficiently  higher  than  one  of  the  other  two  energy  levels, 
Ep  and  Ex. 

The  minima  at  K  =  0  in  the  second  lowest  conduction  mav  have  some 

affect  on  absorption  (and  photogeneration)  at  anal]  wavelengths.  The 

« 

dependence  of  the  energy  of  this  minima,  the  r  minima,  is  given  by 
(Ref  6:3142  ) 

Er-  =  2.90  +  0.36X  +  0.54X2  (86) 

The  pertinent  band  energies  for  each  region  are  indicated  in  Table  I . 
Absorption  Coefficient 

The  saniempirical  absorption  model  of  Hutcby  and  Fudurich  (Ref  6  ) 
is  used  here.  The  constants  introduced  in  this  section  are  also  taken 
from  their  paper. 
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0  .1  .2  .3  .4  .5  .6  .7  .8  .9  1.0 

A1  Content,  X 


Figure  4.  Band  Edges  as  a  function  of  A1  Content,  X, 
in  AlxGa1)(As  (Ref  8;  6:3142) 


The  absorption  coefficient,  a  ,  is  a  function  of  the  wavelength 
of  light,  A,  or  its  corresponding  photon  energy,  E  =  he/ A.  It  is  also 
dependent  on  the  shape  of  the  valence  and  conduction  bands,  in  par¬ 
ticular  the  energies  of  the  minima  in  the  conduction  band.  Therefore, 
a( A)  is  dependent  on  the  A1  content  and  is  in  general  different  in  each 
region . 

The  absorption  coefficient  is  assumed  to  be  the  sum  of  three 
absorption  coefficients  (Ep,  Ex,  and  Er'  respectively). 

a(A)  =  cip  +  ax  +otp'  (87) 

where 

Op  =  Ae  exp  for  E  <  Er 

Ar  •  (E-Ep)1/2  for  E  >  E 

ajj  =  0  for  E  <  Ex 

Ax  •  (E-Ex)2  for  E  >  Ex 

ar'  =  o  for  E  <  Ep' 

Ar,  •  (E-Ep,)1/2  for  E  >  Ep' 

and  where  Ajr  =  1.9  x  102  cm-^ 

Ee  =  1.16  x  10-2  eV 
Ap  =  4.0  x  10^  cm"-*- 
Ax  =  1.0  x  102  cm-'1'  eV~2 
Ap'  =  3.5  x  102  cm"f  eV-^/2 


Relative  Carrier  Effective  Masses 


The  effective  mass  of  a  carrier  is  a  measure  of  its  inertia,  or 
a  measure  of  how  much  energy  is  needed  to  change  its  momentum  The 
effective  mss  of  an  electron  is  determined  by  the  curvature  of  the 
conduction  band  edge,  Ec(k). 

"*  -  2^  <88> 

Ik2 

For  convenience  the  relative  electron  effective  mass  is  defined  as 

Mn  =  (89) 

mo 

where 

o] 

mo  =  electron  rest  mass  (9.11  x  10  kg) 

In  regions  2,  3,  and  4  most  of  the  free  electrons  are  near  the  r 
minima  in  the  conduction  band,  since  it  has  the  lowest  energy.  The 
relative  effective  mass  at  the  r  minima  varies  with  A1  content  according 
to  (Ref  8 ) . 


Mr  =  0.067  +  0.083X  (90) 

In  region  1  most  of  the  free  electrons  are  near  the  X  minima  in 
the  conduction  band.  The  relative  effective  mass  at  the  minima  also 
varies  with  A1  content  (Ref  8). 

%  =  0.85  -0.07X  (91) 

Two  degenerate  valence  bands  are  at  their  common  maximum  at  k  =  0, 
as  previously  mentioned.  The  bands  have  different  curvatures ,  so  the 
holes  in  each  have  different  effective  masses.  The  relative  light-hole 
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effective  mass  (Mih) .  relative  heavy-hole  effective  mass  (M^),  and 
overall  relative  density  of  states  hole  effective  mass  (Mp  =  M^)  are 
independent  of  A1  content  and  are  given  by  (Ref  4:CS-6) 


Mih  =  0.12 


(92) 


Mhh  —  0.62 

3/2  3/2 

U„  =  > 


(93) 

(94) 


Intrinsic  Carrier  Density 

The  intrinsic  carrier  densities  in  regions  2  and  3  are  the  same, 
since  both  regions  have  the  same  A1  content.  The  band  gap  is  Er2 
because  it  is  lower  than  the  other  band  edges.  The  intrinsic  carrier 
density  in  regions  2  and  3  is  (Ref  12:27) 


ni2  =  ni3  =  2  .(2Kg_I)3/2.  (nfe  m*2)3/4  exp  (-E^/2  KBT) 

(95) 


which  can  be  simplified  to 


ni2  =  ni3  =  2.51  x  1019  •  (Mn2  Mp2 )3/4  •  exp  (-Er2/0. 05170  eV) 

•  cm~3 

(96) 

The  diode  current  density,  Jq,  in  Eqn  (73)  is  roughly  proportional  to 
ni2-  ni2  is  a  strong  function  of  Er2.  As  X2  increases,  Er2  increases, 
n^2  decreases,  and  JB  decreases. 

Static  and  Dynamic  Dielectric  Constants 

Harrison  and  Hauser  (Ref  5:298  )  have  developed  a  model  for  the 
static  dielectric  constant,  Ks,  and  dynamic  dielectric  constant,  K^,  for 
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ternary  compounds.  Applying  their  model  to  AlGaAs, 


1  +  2  X  • 


KS)d  - 


„-l 


K^2  +  a-X) 


vP  ,-l 
s,d 


1  -  x  •  ^.d-1 


(i-x)  j&jli 

KS)d+2 


where 


Ks  =  static  dielectric  constant  of  AlGaAs 


K(j  =  dynamic 
Kg  =  static 
=  dynamic 
Kg  =  static 
K§  =  dynamic 


"  AlGaAs 


"  AlAs 


"  GaAs 


The  following  dielectric  constants  are  used  for  AlAs  and  GaAs  (Ref 
11:75): 


Kg  =  10.9 


Kg  =  13.2 


=  8.5 
K§  =  10.9 


Using  these  constants  in  Eqn  (97)  gives  the  dielectric  constants  for 
AlGaAs  as  a  function  of  A1  content. 


„  _  2.60  -  0.70  X 

cyr?  n  7P5R  v 


x,  _  2.54  -  0.106  X 
M  ~  "0.233  -  0.053“X 
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and  Ionization  Energies 


Dopants 

The  ionization  energy,  Ei,  of  a  dopant  in  a  particular  semi¬ 
conductor  determines  the  fraction  of  dopant  atoms  ionized.  Be 
(Beryllium)  is  used  as  the  p  dopant  in  regions  1  and  2,  and  Te 
(Tellurium)  is  used  as  the  n  dopant  in  regions  3,  4,  and  5.  Be  has 
an  ionization  energy  of  0.030  eV  in  GaAs  (Ref  7).  Te  has  an  ioniza¬ 
tion  energy  of  0.003  eV  in  GaAs  (Ref  13).  The  ionization  energy  is 
roughly  proportional  to  the  conductivity  effective  mass  (of  conduc¬ 
tion  electrons  for  donors,  and  of  holes  for  acceptors)  and  inversely 
proportional  to  the  square  of  the  static  dielectric  constant  (Ref  12: 
29) .  The  conductivity  effective  mass  and  density  of  states  effective 
mass  are  approximately  equal  for  III  -  V  semiconductors.  The  ioniza¬ 
tion  energy  of  Te  as  a  donor  in  AlGaAs  is  given  by 


where  G  indicates  GaAs.  More  simply,  for  regions  3,  4,  and  5 


St'  (-003  eV> 


(101) 


Since  the  hole  effective  mass  is  independent  of  A1  content  and  even 
simpler  result  applies  to  the  ionization  energy  of  Bg  as  an  acceptor 
in  region  1  and  2. 


Ei 


(.030  eV) 


(102) 


Doping  Levels 

As  already  mentioned,  the  solar  cell  has  been  designed  to  con¬ 
strain  minority  carriers  to  flew  toward  the  junction.  However,  it  is 
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also  desirable  to  allow  majority  carriers  to  flow  easily  away  from 
the  junction.  If  the  potential  energy  of  electrons  in  region  4  is 
higher  than  in  region  3,  the  electrons  will  not  move  easily  frcm 
region  3  to  region  4.  This  problem  would  show  up  as  an  internal 
resistance.  The  electron  thermal  equilibrium  concentration  in  region 
4  is  required  to  be  greater  than  or  equal  to  that  in  region  3. 

N40  >  N30  (103) 

Since  regions  3  and  4  are  in  contact  seme  electrons  from  region  4 
will  cross  the  boundary,  causing  a  net  negative  charge  near  the 
boundary  in  region  3  and  a  net  positive  charge  in  region  4.  Conse¬ 
quently,  an  electric  field  exists  across  the  boundary,  which  allows 
easy  flow  of  electrons  frcm  region  3  to  region  4. 

In  order  to  derive  the  requirement  on  the  total  concentration 
of  donors  (ionized  and  unionized)  in  regions  3  and  4,  the  following 
definitions  are  made: 
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Ej_  =  -  Eq  =  donor  ionization  energy 

N*  =  -  Nc  e_Ei/fkBT 
g  <- 


The  concentration  of  unionized  donors  is  given  by  the  Fermi-Dirac 
distribution. 


So, 


N°  -  Nfc 
D  D 


1  +  -  e(ED-EF)/kBT 

g 


^  “  i+e (Ec~Ep) /kgT  m  e (Ed-Ec) /kg’! 

Nc  *  g 


(104) 


(105) 


But, 


Nq  =  Nc  e" (Ec4-Ef) /kgT 


(106) 


Therefore , 


Nd  *d  [  1  ~  l  +  N*/No 


Nf 


Nn  N* 
Nq+N* 


(107) 

(108) 


The  thermal  equilibrium  concentration  can  also  be  written  as 


Nq  = 
No  = 

No  = 
Nq  =  (No 
<No 


t  n1 

Po  +  ND  =  t 


+  Ni 


2  Nn  N* 
ni  +  -Jd - 

Nq  ^ 


n2  (Nq  +  N*)  +  Ng  N*  Nq 
Nq(N0+N*) 

+  N*)  =  n?  (No  +  N*)  +  n£  N*  Nq 
-  n?)  (N0  +  N*)  =  N°  N*  N0 


(109) 

(110) 

(111) 

(112) 

(113) 
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(114) 


For  high  doping  levels  N§  >>  n?  so 


Nq  (Nd  +  N*)  =  N*  Nq 


Solving  for  NQ 


„  -  N* 
Nq - 2 


O 

Eqn  (115)  applies  to  regions  3,  4,  and  5.  is  chosen  to  be 

1  x  1018  cm-3,  which  implies  N4q  =  4.139  x  10^7  cnf3.  Eqn  (103)  is 
easily  satisfied  for  <  1  x  1018  cm'3  and  X3  <  .3  as  demonstrated 

in  the  following  Table  for  =  1  x  lO^8  cm-3 


X 

N30  (cm“3) 

0 

3.538  x  1017 

.1 

3.758  x  1017 

.2 

3.959  x  1017 

.3 

4.139  x  1017 

A  requirement  similar  to  Eqn  (103)  applies  regions  4  and  5. 

N50  >  N40  (116) 

Ng  is  chosen  to  be  1.5  x  103-®  an-3,  which  implies  N50  =  4.508xl03-7cm-3, 
so  that  this  condition  is  satisfied. 

The  thermal  equilibrium  hole  concentrations  in  regions  1  and  2 
are  determined  by  an  equation  similar  to  Eqn  (115) . 


p*  / 4n£  ' 

Po  =  f-Vl+-p£  -1 

where 
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=  total  concentration  of  acceptor  atoms 
p*  -  ^  e  -Ei/kBT 

g 


g  =  4  (two  degenerate  bands  x  2  spins) 

Nn  =  2  (£  =  density  of  states  in  the  valence  band 

h2 

The  doping  requirement  for  regions  1  and  2  is 


P2o  >  Pio  (H8) 

is  chosen  to  be  1.5  x  lO^8  cm-8,  which  implies  P^0  = 

7.095  x  1017  cm-8.  Eqn  (118)  is  then  easily  satisfied  for  < 

1  x  1018  an-8.  In  all  computer  runs  P2q  and  N30  were  printed  out  to 
insure  that  Eqns  (103)  and  (118)  were  satisfied. 


Minority  Mobilities  and  Diffusivities 

The  model  of  Sutherland  and  Hauser  (Ref  11:62-64,  76)  is  used 
in  this  section  to  determine  the  minority  carriers  mobilities.  The 
model  is  based  on  the  assumption  that  mobility  is  primarily  deter¬ 
mined  by  a  polar  optical  phonon  scattering  process. 

In  region  2,  virtually  all  conduction  electrons  are  in  the 
direct  valley  (r) .  This  is  also  true  of  GaAs.  The  electron  mobility 
in  region  2  is  then  related  to  the  electron  mobility  in  GaAs  (G 
indicates  GaAs) . 
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Diffusion 


Diffusion  lengths  are  very  fabrication  dependent.  Even  when 
using  the  same  technique  and  same  equipment  experimenters  obtain 
samples  with  widely  varying  diffusion  lengths.  For  this  reason 
a  conservative  value  of  2y  will  be  used  in  all  regions. 

Surface  Recombination  Velocities 

Surface  recombination  velocities  are  also  fabrication  dependent. 


Therefore,  they  are  left  as  variables  in  the  program.  A  typical 
value  is  10 4  cm/ sec.  For  convenience  all  three  surface  recombina¬ 
tion  velocities  will  be  varied  together  (S  =  SQ  =  Sp  =  S3) . 


V.  Computer  Model  and  Results 


A  description  of  the  program  used  to  model  the  AlGaAs  solar  cell 
is  described  in  this  chapter.  The  results  of  varying  seme  of  the 
parameters  are  presented  and  discussed.  It  should  be  noted  that 
this  chapter  does  not  include  the  entire  results  of  this  report. 

In  particular,  the  theoretical  results  discussed  in  Chapter  II  are 
fairly  general  and  include  important  statements  on  multi-layered 
solar  cells. 

Description  of  Program 

The  program  used  to  model  the  cell  is  based  on  the  results  of 
the  previous  three  chapters.  The  inputs  to  the  program  are  X 
(=  X2  -  X3) ,  Z2,  Z3,  and  S.  Other  parameters  were  left  as  variables 
in  the  program  that  can  easily  be  changed  or  made  inputs  to  the 
program.  However,  in  the  present  study  the  four  parameters  mentioned 
above  are  of  primary  interest.  The  other  variables  were  held  at  the 
following  values 

Z^  =  . 3y 
Z4  =  l.Ou 

N°2  =  1  x  1018  cm-8 

N°  =  1  x  lO*8  cm-8 
D3 

Li  =  L2  =  L3  =  L4  =  2u 

The  program  used  the  inputs,  the  above  variables,  values  for  parameters 
specified  in  Chapter  IV,  and  the  Eqns  in  Chapter  IV  to  calculate  all 
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the  parameters  needed  in  Chapters  II  and  III.  Table  I  in  Chapter 
IV  indicates  many  of  the  parameters  specified  or  given  in  terms  of 
other  parameters .  These  were  then  used  in  Eqns  (61)  through  (64) 
in  Chapter  II  to  compute  the  spectral  current  densities.  The  inte¬ 
gration  indicated  in  Eqn  (66)  was  done  numerically  using  wavelength 
increments  of  AX  =  .Olp.  Jg  c  was  obtained  as  in  Eqn  (67).  JD 
was  determined  using  Eqns  (70)  through  (73) .  Eqn  (74)  was  used  to 
determine  JT  as  a  function  of  V,  with  increments  of  V  =  .001  volts. 

At  each  V  the  efficiency,  n,  was  calculated  from  Eqn  (78).  These  were 
compared  and  the  maximum  n,  the  efficiency  of  the  cell,  found — 
along  with  its  corresponding  JT  -  V  pair.  The  open  circuit  voltage 
Vqc  was  obtained  by  interpolating  between  the  points  at  which  JT  was 
just  above  and  just  below  2ero. 

Results 

At  S  =  10^  cm/sec  (typical  value)  the  A1  content  in  regions  2 
and  3  (X  =  X2  =  X3)  was  varied.  At  X  =  0,  .05,  .01,  and  .02,  Z2  and 
Z3  were  varied  and  the  optimum  (for  maximum  efficiency)  combination 
of  Z2  and  Z3  found.  At  each  of  these  A1  contents  the  optimum  values 
of  Z2  and  Z3  were  found  to  be  Z2  =  0.50p  ±  .Olp  and  Z3  =  2.4p  ±  .lp. 

X  was  then  varied  in  small  increments  from  0  to  .25  and  some  of 
the  characteristic  of  the  cell  noted.  These  characteristics  are 
presented  in  Table  II.  The  efficiency,  n>  was  found  to  be  relatively 
independent  of  X  for  0  <  X  <  ,15.  The  optimum  value  of  X  was  X  =  .05 
where  n  =  17.71%.  However,  this  is  only  .23%  higher  than  for  X  =  0 
(plain  GaAs)  where  n  =  17.58%.  n  varied  in  a  complicated  manner 
with  X  over  the  range  0  <  X  <  .12  achieving  several  maxima,  n  de¬ 
creased  with  increasing  X  for  X  >  .12. 
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17.269  5.363  22.632 


rn 

The  short  circuit  current,  c  ,  decreased  steadily  with 
increasing  X,  while  the  open  circuit  voltage,  Voc<,  increased 
steadily  with  increasing  X.  This  is  as  anticipated.  The  band  gap 
increases  with  X  and  leads  to  smaller  absorption  coefficients 
{especially  at  large  \) ,  and  a  smaller  intrinsic  carrier  density  n^ 

(see  Table  II) .  The  smaller  absorption  coefficients  lead  to  less 
photogeneration  particularly  at  large  A.  This  is  evident  by  comparing 
Figures  5,  7,  and  9,  where  the  spectral  current  densities  for  X  =  0, 

.05,  and  .25  are  presented.  The  JT  versus  V  plots  for  these  values 
of  X  are  given  in  Figures  6,  8,  and  10. 

The  A1  content  in  regions  2  and  3  was  also  varied  at  S  =  0  and 
S  =  10?  an/ sec.  For  S  =  0  {no  surface  recombination)  the  optimum 
values  of  Z2  and  Z3  were  found  to  be  Z2  =  0.48u  and  Z3  =  1.8u 
independent  of  Al  content.  The  results  for  S  =  0  are  given  in  Table 
III.  They  are  similar  to  those  for  S  =  10^.  Again  the  efficiency 
varied  in  a  complicated  manner  with  X,  achieving  several  local  maxima 
in  the  range  0  <  X  <  . 12  but  decreasing  for  X  >  . 12 .  The  best  ef¬ 
ficiency,  18.01%,  is  again  just  slightly  higher  than  at  X  =  0  (17.87%). 
Efficiencies  were  slightly  higher  for  S  =  0  than  for  S  -  10^  cm/ sec 
because  more  of  the  photogenerated  minority  carriers  reach  the  junction. 

The  results  for  S  =  10^  cm/sec  are  given  in  Table  IV.  The  ef¬ 
ficiency  was  found  to  increase  steadily  with  decreasing  Z2  and  in¬ 
creasing  Z3  for  all  X  used.  A  small  Z2  and  large  Z3  are  otpimum  for 
very  large  S,  because  most  of  the  photogenerated  minority  carriers 
will  then  be  produced  in  region  3  near  the  junction,  where  the  only 
surface  recombination  that  can  affect  them  is  far  away  at  y3. 
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WAVELENGTH  (/*) 


Fiqure  5.  Photogenerated  Spectral  Current  Densities 
for  Sa]o^  cm/sec,  Xa0,  Z2*0  •  50j»  ,z  ^*2 . 4^ 
(optimum  values  Z£  and  Z3) 


49 


VOLTAGE(volts) 

0T  versus  V  for  S«10*cm/sec ,  X*.05, 
22*0*50 ja,  Zj«2.4 j*  (optimum  Z£  and  z3) 


s 

2 

3 
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Table  III 


Solar  Cell  Characteristics  for  Different  Al  Content  in  Regions  2  and  3 
(X  =  X2  =  X3)  at  S  =  0,  and  Optimum  Z2  and  Z3  (Z2=0.48u,  Z3=1.8y) 


-rT 

'■’s.c. 

Vo.c. 

n 

X 

(mA/cxn^) 

(Volts) 

(%> 

0  33.265  .9266  17.87 


.01 

32.865 

.9389 

17.93 

.02 

32.415 

.9511 

17.96 

.03 

31.956 

.9633 

17.98 

.04 

31.476 

.9755 

17.97 

.05 

31.073 

.9878 

18.01 

.06 

30.629 

1.0000 

18.01 

.07 

30.110 

1.0120 

17.96 

.08 

29.768 

1.0243 

18.01 

.09 

29.337 

1.0365 

18.00 

.10 

28.848 

1.0485 

17.94 

.11 

28.504 

1.0608 

17.97 

.12 

28.080 

1.0730 

17.94 

.13 

27.658 

1.0851 

17.90 

.14 

27.274 

1.0973 

17.89 

.16 

26.48] 

1.1215 

17.82 

.18 

25.691 

1.1457 

17.72 

.22 

24.048 

1.1938 

17.40 

.25 

23.127 

1.2304 

17.33 
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Table  IV 


Efficiency  as  a  Function  of  Al  Content  in  Regions  2  and  3 
(X=X2=X3)  at  S=107cnv'sec/  and  22=0. 20y  and  Z3=10u 


n 


X 

(%) 

0 

12.69 

.01 

12.71 

.02 

12.70 

.03 

12.68 

.04 

12.65 

.05 

12.65 

.06 

12.66 

.07 

12.55 

.08 

12.57 

.09 

12.53 

.10 

12.46 

.12 

12.42 

.15 

12.28 

.20 

11.98 

.25 

11.69 
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n  increased  only  slightly  (.01%)  for  Z3  above  lOu  so  this  value  was 
chosen  for  comparison.  Since  it  is  difficult  to  make  layers  smaller 
than  . 2u  this  value  was  used  for  Z2-  The  efficiencies  were  much 
less  than  for  S  =  0  or  S  =  10^  cm/ sec.  Again  little  was  gained  by 
increasing  Al  content.  In  fact,  for  X  >  .03  n  was  less  than  for 
X  =  0.  The  optimum  value  of  X  =  .01  gave  only  a  .02%  increase  in  n 
above  that  for  X  =  0  (12.71%  versus  12.69%). 
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VI. 


and  Conclusions 


The  multi-layered  cell  indicated  in  Figure  1  (a  solar  cell  with 
a  window,  a  single  hcmo junction,  and  a  back-barrier  on  a  base)  has 
been  analyzed  theoretically  under  fairly  general  conditions.  The 
theoretical  results  are  applicable  to  similarly  structured  cells 
made  from  other  materials.  The  carrier  transport  equations  were  solved 
to  obtain  the  excess  minority  carrier  densities  and  the  short-circuit 
current  density.  The  back-barrier  was  proven  theoretically  to  in¬ 
crease  the  short  circuit  current  and  therefore  the  efficiency  of  the 
solar  cell.  The  back-barrier  is  most  useful  when  low  surface  recombina¬ 
tion  velocities  can  be  obtained,  and  when  a  large  portion  of  the  photo¬ 
generation  takes  place  in  the  region  in  front  of  it.  The  program 
already  developed  could  be  easily  modified  to  model  a  cell  without  a 
back-barrier.  The  increase  in  efficiency  due  to  the  back-barrier 
could  then  be  obtained  for  various  conditions. 

The  properties  of  AlxGa^_xAs  have  not  been  as  thoroughly  in¬ 
vestigated  as  those  of  other  materials.  This  problem  is  compounded 
because  data  on  these  properties  is  needed  as  a  function  of  A1  con¬ 
tent.  In  this  report  many  of  the  properties  of  AlxGa-^_xAs  have  been 
modeled  by  drawing  on  the  work  of  others.  Although  much  research  on 
AlGaAs  is  still  needed,  Chapter  IV  includes  a  significant  up-to-date 
ccmpilation  of  information  on  and  models  of  the  properties  of  this 
material . 

The  AlGaAs  solar  cell  was  modeled  on  a  computer  and  sane  of  the 
design  parameters  varied  in  order  to  examine  their  effects.  The 
optimum  widths  of  regions  2  and  3  were  found  to  be  virtually  inde¬ 
pendent  of  A1  content,  X,  for  0  <  X  <  .23,  but  dependent  on  the  surface 
recombination  velocity. 
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Only  small  gains  in  efficiency  (less  than  0.3%)  were  obtained 
by  varying  the  A1  content  in  regions  2  and  3.  However,  these  results 
were  based  on  the  assumption  that  the  diffusion  lengths  were  inde¬ 
pendent  of  Al  content.  AlGaAs  can  be  expected  to  have  a  higher 
diffusion  length  than  GaAs  when  made  from  liquid  phase  epitaxial 
growth,  since  the  Al  in  the  melt  acts  as  a  gatherer  of  impurities. 
This  increase  in  diffusion  length  would  result  in  a  slightly  higher 
short  circuit  current  and  a  significant  improvement  in  voltage  per¬ 
formance  and  fill-factor  (see  Eqn  (73) ) .  Also  the  cell  analyzed 
can  be  used  in  a  concentrator  application  wherein  the  voltage  operat¬ 
ing  point  would  be  such  that  the  thermal  diffusion  dominates  the  dark 
current-voltage  characteristic  resulting  in  a  significant  increase  in 
fill-factor. 

Since  Air  Force  missions  in  space  are  limited  to  earth  orbit, 
modifications  of  this  cell's  design  will  be  required  to  optimize  its 
performance  in  radiation  belts  (electrons  and  protons)  for  flat  panel 
applications.  However,  very  little  data  exists  on  radiation  damage 
effects  to  AlGaAs  material. 
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